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Abstract 

We propose a non-perturbative method for defining the higher dimensional opera¬ 
tors which appear in the Heavy Quark Effective Theory (HQET), such that their matrix 
elements are free of renormalon singularities, and diverge at most logarithmically with 
the ultra-violet cut-off. Matrix elements of these operators can be computed numer¬ 
ically in lattice simulations of the HQET. We illustrate our procedures by presenting 
physical definitions of the binding energy (A) and of the kinetic energy (-Ai/2mQ) of 
the heavy quark in a hadron. This allows us to define a “subtracted pole mass”, whose 
inverse can be used as the expansion parameter in applications of the HQET. 
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1 Introduction 


The Heavy Quark Effective Theory (HQET) has developed during the last few years 
into a very useful tool for the study of strong interaction effects in heavy quark physics 
i-i (for a comprehensive review see ref. [^ ). In this approach physical quantities are 
studied systematically as series in inverse powers of the mass of the heavy quark [|. In 
particular, local composite operators of QCD, whose matrix elements contain the long¬ 
distance gluonic effects in physical processes, are expanded in terms of the operators 
of the HQET 

where mg is the mass of the heavy quark and // is the renormalisation scale used 
in defining the renormalised operators Eor the moment we assume that 

the operators are renormalised in some scheme based on the dimensional 

regularisation of ultra-violet divergences, such as the MS scheme. The QCD opera¬ 
tors generally also depend on some renormalisation scale, which we take to be 

different from //, and we do not exhibit this dependence explicitly. Renormalon singu¬ 
larities lead to ambiguities in the Wilson coefficient functions and in the matrix 

elements of the operators Although these ambiguities cancel (see ref.Q and 

references therein, and refs. i-H), the presence of renormalons requires an alterna¬ 
tive definition of the renormalised operators if the HQET is to be applicable 

beyond leading order in the heavy quark mass, and in particular if the coefficient func¬ 
tions are to be calculable using perturbation theory. The renormalon singularities and 
the corresponding ambiguities which we are considering here are those induced by the 
introduction of the expansion in eq.(j^. Of course the matrix elements of the QCD 
operator will themselves contain non-perturbative long-distance effects, but, as 

usual, these effects do not appear in the coefficient functions, but only in the matrix 
elements of the operators of the HQET. 

The use of a hard (dimensionful) ultra-violet cut-off A in the effective theory (in¬ 
stead of dimensional regularisation) leads to matrix elements and coefficient functions 
which are free of renormalon ambiguities |^] . However in this case the matrix elements 
in the effective theory diverge as powers of the cut-off A. The subtraction of these 
power divergences cannot be performed using perturbation theory; such a subtraction 
reintroduces renormalon ambiguities in the matrix elements and coefficient functions. 

This point will be discussed in detail below. 

^In B D decays the expansion may be performed simultaneously in the inverse powers of the masses 
of both heavy quarks. 


qQCD ^ 

n,Q. 
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In this paper, in contrast to previous approaches, we propose a non-perturbative 
definition of the renormalised operators such that they are free of both renor- 

malon ambiguities and power divergences. This allows us to present “physical” defini¬ 
tions of parameters such as the binding energy (A) or the kinetic energy (—Ai/2mQ) of 
the heavy quark in a heavy hadron. By “physical” we mean that they do not depend 
on any regularisation or renormalisation scale, nor on the method used to regulate the 
ultra-violet divergences (although they do depend on the renormalisation prescription 
used to define them). 

An immediate, and related, question concerns the choice of the mass parameter tuq 
used in the expansion. The pole mass is ambiguous due to the presence of infra-red 
renormalon singularities ||^, |^, and so a different definition of niq is required. We 
introduce the definition of a “subtracted pole mass”, mg, from which the renormalon 
ambiguities have been subtracted non-perturbatively. mg has the attractive property 
that it can be determined from the mass of a hadron H containing the heavy quark 
Q, by calculations performed entirely within the HQET. For each such hadron H we 
define a parameter A = mH — and A — > constant as mg —> oo (for simplicity of 
notation we suppress the label H on K). A remains finite when the ultraviolet cut-off is 
taken to infinity, and can be computed numerically in lattice simulations of the HQET. 

We discuss the question of the definition and evaluation of A to leading order in section 
^ and to 0{l/mQ) in section]^ 

As a further illustration of our general procedure we study in section ^ the matrix 
elements which determine the kinetic energy of the heavy quark in a hadron ^ 

^ {H\h&h\H) 

= - 

2mH 

where h represents the heavy quark field in the HQET, and H is a hadron (of mass mn) 
containing a heavy quark (we suppress the label H on Ai). The Ai’s are important 
ingredients in the study of the spectroscopy and inclusive decays of heavy hadrons. 

The chromomagnetic operator haijG^^h, where Gq are the spatial components of the 
gluon field strength tensor, appears at the same order of the heavy quark expansion 
as the kinetic energy operator. However the spin structure of the chromomagnetic 
operator ensures that its matrix elements are free of renormalon singularities. 

In ref.|]T^, together with L.Maiani, we pointed out that in the lattice formulation 
of the HQET the evaluation of higher order terms in l/rug involves the appearance 

^Here we consider the effective theory for a heavy quark at rest. The generalisation to an arbitrary 
four-velocity is straightforward, although there are considerable subtleties in formulating the effective theory 
at non-zero velocity in Euclidean space BID- 
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of ultra-violet singularities which diverge as inverse powers of the lattice spacing a. 
This is due to the fact that the higher order terms generally involve matrix elements of 


higher dimensional operators which can mix with lower dimensional ones (e.g. hD^h 
can mix with l/ahDi^h and 1/a^hh). In ref.]^ it was further stressed that these 
divergences must be subtracted non-perturbatively, since factors such as 



(3) 


which do not appear in perturbation theory, give non-vanishing contributions as a —> 0. 
Renormalons represent an explicit example of non-perturbative effects of this kind. In 
this paper we explain how the subtraction of power divergences and of renormalon 
singularities, can be performed non-perturbatively, illustrating our ideas with the eval¬ 
uation of A and Ai. 

Much of the discussion in this paper is presented within the framework of lattice 
field theory, however all the theoretical questions addressed below have to be faced 
with any ultra-violet regularisation scheme ||^ . Moreover, our proposed definitions of 
subtracted higher dimensional operators in general, and of A and Ai in particular, are 
in fact independent of the regularisation. The primary aim of this paper is to provide 
an understanding of these theoretical issues. In addition however, lattice simulations 
provide the opportunity for the evaluation of the parameters of the HQET, such as A 
and Ai, from first principles. In the following sections we also explain how this can be 
done in principle, and in ref.|]l^ we present some results from an exploratory numerical 
study (some preliminary numerical results, and a summary of the ideas of this paper 
have been presented in ref.[^]). However, instead of using lattice simulations, it is also 
possible to use other non-perturbative methods, such as QCD sum rules, to compute the 
matrix elements of the operators defined by our prescription. These matrix elements 
are free of renormalons and power divergences. 

Throughout this paper we use the following notation for ultra-violet cut-offs and 
renormalisation scales. We denote by M the scale used to define the renormalised 
operators in QCD and by /i that to define the operators in the HQET (the 

Thus ^ -C mQ <C M. It will also be convenient to consider bare operators 
in the HQET, and we denote by A the corresponding ultra-violet cut-off. Where the 
discussion is particular to the lattice formulation of the HQET we replace A by a~^. 

The plan of the remainder of this paper is as follows. In the next section we 
discuss the definition of A, and explain how it can be computed in lattice simulations. 
We explain our general procedure for the definition of subtracted higher dimensional 
operators in section using the matrix elements Ai as examples. Section Q contains 
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a discussion of “matching”, i.e. the determination of the coefficient functions Cn,a 
corresponding to our definition of the operators. Finally in section ^ we present our 
conclusions. 


2 The Binding Energy - A 

In this section we discuss the definition of the heavy quark mass tuq, from which 
the ambiguities from the leading renormalon singularities have been subtracted non- 
perturbatively. mg can thus be used as the expansion parameter in the HQET. We 
start however by reviewing very briefly some of the relevant points concerning the 
renormalon singularities in the propagator of the heavy quark. 


2.1 Renormalons in the Heavy Quark Propagator 


Our discussion follows closely the presentation by Beneke and Braun Q, in which 
these authors study the heavy quark propagator in the large N^ limit, where Nf is the 
number of light quark flavours. We refer the reader to ref. pi] for more details, and also 


to ref. [1^1 in which the discussion is extended to include the terms of ©(l/mg) in the 
inverse propagator. To leading non-trivial order in 1/Nf, the renormalon singularities 
are obtained by summing over an arbitrary number of light quark loops in the gluon 
propagator. The Borel transformed gluon propagator (in the Landau gauge) is written 


as 









- k'^g, 


flU 


(_^2)2+n 


(4) 


where Pq is the coefficient of a^/dTr in the first term of the /3-function (/3o = — (H — 
2/3Nf)), and A and B are colour labels. // is the renormalisation scale and C is a 
scheme dependent constant. is the contribution to the gluon propagator from 

the diagram with n quark loops. Of course in leading order in Nf it is only the term 
2/3Nf which appears, however it is assumed that the replacement of 2/3Nf by Pq 
is a consistent one for identifying the singularities. This assumption is based on the 
intuition, that it is the infra-red behaviour of the running coupling constant which is 
(at least partially) responsible for the singularities. 

Feynman diagrams evaluated with the Borel transformed propagator (Q) may have 
poles at positive values of u, rendering the inverse Borel transform (which requires an 
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integral over positive values of u) ambiguous |^. This is particularly true in Operator 
Product Expansions, where such singularities in the Wilson coefficient functions are 
cancelled by those in the matrix elements of higher dimensional, or higher twist, op¬ 
erators 1^, 1^. The expansion in inverse powers of the mass of the heavy quark is an 
interesting example of this phenomenon. We start by considering the quark propagator. 
The inverse quark propagator in QCD can be written in the form 

S~^ {p,m) =jjl — m — Y,{p,m) (5) 

where 

S(p, m) = mSi(p^, m) + (^ — m)Y, 2 {p‘^,m) (6) 

and m is the bare mass. We now write 


p = niQV + k 


(7) 


where ruQ is some well defined choice of the heavy quark mass, to be specified later. 
In all the explicit examples given in sections |2.2| and |3| we will take u = (1,0). It is also 
convenient to define the quark propagator sandwiched between projection operators 


1 +y 
2 


Sp{k,mQ) 


i+i 

2 


S{p, m) 


i+i 

2 


( 8 ) 


Then, [§, 


Sp^{k,mQ) 


mQ-mpo\e{mQ, p)-C{mQ/p) S^g{v-k, p)+0 


/ {v ■ k)'^ k‘j_ 1 \ 
V mg ' mg' Nj) 


(9) 


where 5efT is the quark propagator in the HQET (whose action is given by hiv ■ Dh] 
h represents the field of the heavy quark), p is the renormalisation scale, and k^ = k"^ — 
{vkf. s:i is linearly divergent in perturbation theory. In dimensional regularisation 
the Borel transforms of both rUpoie and S~^ have renormalon singularities at tt = 1/2, 
the infra-red renormalon in mpoie cancels the ultra-violet renormalon in S~^ on the 
r.h.s. of eq. (® 0- The connection between power divergences and renormalons can be 
seen by noting that if S~q is linearly divergent in one-loop perturbation theory (which 
corresponds to using (^) with u = 0 as the gluon propagator), then its Borel transform 
is logarithmically divergent at u = 1/2, corresponding to a pole singularity at this 
point. In particular S'/g^(0,/x) is not equal to zero, which is a signal of the presence 
of ultra-violet renormalon singularities at u = 1/2. A more detailed account of the 
correspondence between power divergences and renormalons is given in ref. HI- 

^When the calculations are extended beyond the leading order in 1/lV/, the poles become replaced by 
branch points of cut singularities. 
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For the discussion below, it is convenient to re-express eq. (I) in terms of the bare 
propagator in the effective theory with a hard cut-off A, 


Sp^{k,mQ) 


mg - mpoie(mQ, A) - C{mQ/A) ■ k,A) + 0 


/ {v ■ k^ 1 \ 

V ^ rng' Nj) ' 
( 10 ) 


Now mpoie(mQ, A) and S~^{v ■ k,A) both diverge linearly with the cut-off, but their 
Borel transforms have no poles at u = 1/2. Thus we have replaced the problem 
of the ambiguities associated with renormalon singularities with that of determining 
the exponentially small terms when power divergences are present (see eq.(|^ and 
the corresponding discussion). By adding a residual mass counterterm {6m hh) to the 
action of the HQET, it is possible to impose the non-perturbative condition *S'/^^(0, A) = 
0, thus removing the power divergences. In the following subsection, we impose instead 
an equivalent condition on the static (i.e. F = 0) heavy quark propagator in Euclidean 
configuration space {S{^g{x-,y), where i and j are colour labels) in the Landau gauge: 


^lii(l‘S'eff(^;0)|) = 0 (11) 

where t = and is the trace over the colour components of the propagator. Since 
the operator we are subtracting {hh) is conserved, either of these conditions is sufficient 
to determine its coefficient 5m fully. Matching the effective theory onto full QCD now 
implies that mpoie on the r.h.s. of eq.([IoD is replaced by a “subtracted pole mass”, 
ttiq, from which the power divergences (and indeed all the dependence on A) have 
been subtracted, mg is therefore a natural choice for the expansion parameter mg of 
the HQET. The implementation of the above subtraction requires the non-perturbative 
evaluation of the quark propagator in the HQET because of the arguments given in the 
discussion of eq.(P). Lattice simulations provide the possibility for such an evaluation, 
and in the following subsection we discuss an explicit definition of mg based on eq. (|TT|) , 
and of the corresponding A parameter, A = mn — '^g: where H is a hadron containing 
one heavy quark Q. In spite of the fact that the entire discussion of the following 
subsection is in the framework of lattice field theory, the value of A is determined 
by its definition based on the condition and is independent of the method of 

regularisation. 


2.2 Definition of A 

In this subsection we present our definition of A explicitly, and also discuss how it might 
be evaluated using the lattice formulation of the HQET. The discussion is presented in 
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Euclidean space, and we take as the action of the HQET in the static (i.e. v = 0) case 


jCefT — hDii^h 


( 12 ) 


where h represents the field of the heavy qnark. Consider the correlation function 

C(t) = ^(0|Jr(f,t) 4(0,0)10) (13) 

X 

where 4 Jr are interpolating operators which can create or annihilate a meson 
state in the HQET. For example, we may take J = hTq where q represents the field 
of the light quark and E is one of the Dirac matrices ^ For sufficiently large times, so 
that only the gronnd state contributes significantly to the correlation fnnction, 

C{t) Z^exp(—ft) (14) 


where the constant Z is independent of the time. The exponent £ is equivalent to the 
definition of A proposed in ref.|l8| 


_ -^4 (0| hrq \M) 
{0\hrq\M) 


(15) 


However £, and hence A defined through eq.([T^), is not a physical qnantity, since it 
diverges linearly as o ^ 0 (as can be demonstrated in one-loop perturbation theory). 
If the matrix elements in eq.(^) are defined in the MS scheme, then the correspond¬ 
ing definition of A contains an ambiguity of 0{Aqcd)- The linear divergence in £ 
appears due to the mixing of the operator JiD^h in the effective action, with the lower 
dimensional operator hh. It is possible to subtract this divergence non-pertnrbatively 
by adding a residual mass term 6m hh to the action, where 6m is determined by im¬ 
posing a suitable renormalisation condition. We now state explicitly onr proposal for 
a physical definition of A, and for its non-perturbative evaluation: 


i) Evalnate the heavy qnark propagator in the theory defined by the action (12), 
in some fixed gauge (the Landau gauge say). In practice this will be done using 
lattice simulations. The heavy qnark propagator in a smooth gange, such as the 
Landau gauge, is of the form (for t > 0) 

{x, f; 0,0) = (0 I W{x)h\Q) | 0) = 6'^{x) 6^^A{t) exp(—At) (16) 


where the ultra-violet divergences associated with the quark mass are contained 
in the exponent A. i and j are colour labels, and unless specifically reqnired they 
will be suppressed below. A{t) satisfies the condition 

lim ^ln{|A(f)|} = 0 (17) 

_ t^oo dt 

^The generalisation to heavy baryons is straightforward. 
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In the explicit examples below we will use the following lattice covariant derivative 

Difix, t) = ^ (/(x, t) - Ul {x, t - a) fix, t - a)j (18) 

where t)} are the link variables. Other lattice definitions of Di are also 

acceptable. 

ii) Add a residual mass term to the effective action 

^eff = ?- (hDih + 5mhh) (19) 

1 + oma 

where 5m will be specified in iii) below. With the new action, the heavy quark 
propagator (now denoted by S'yf) is given by 

S'y^ix, t] 0,0) = Shix, t; 0,0) exp(— ln(l + 6m a)t/a) (20) 


This result demonstrates that the divergences associated with mass renormalisa¬ 
tion do indeed exponentiate. 

iii) We fix the residual mass counterterm 6m by the condition 


ln(l + (lmo) 1, ( \Sl^ix,t + a;6,\ 

— V = - = iim — In I 


t—>oo a 


|S';^*(x,t;0,0)| 


( 21 ) 


u can be computed in numerical simulations by studying the logarithm on the 
right hand side of eq.(^) as a function of t. 

iv) The “physical” definition of A is 


K = £ — V = E + — ln(l + 6m a) 
a 


( 22 ) 


A is finite in the limit a —> 0 and is free of renormalon ambiguities. The corre- 

(23) 


spending “subtracted” pole mass mg is defined by 


mq = mn — A . 


mq is independent of a, and contains no renormalon ambiguities. Thus from the 
computed value of mq one can determine m-j^i^), or any other short-distance 
definition of the quark mass (up to uncertainties which are now of Oi\/mq)), 
using perturbation theory. This will be discussed in section 
The counterterm v defined above is gauge-invariant, in spite of the fact that it is 
calculated from the heavy quark propagator in a fixed gauge. The argument goes 


®The normalisation factor 1/(1 -|- Sma) is introduced for convenience, as will become apparent below. 







as follows. The linear divergence is eliminated from any correlation function, i.e. 
for any external state, by subtracting from the action ( 12 ) a term proportional to 
the gauge-invariant operator hh. Since in this way one eliminates all divergences 
both for quark and hadron external states, the coefficient of the mixing has to 
be gauge-invariant. This must be true also for the finite non-perturbative term 
which accompanies the linear divergence. The gauge-invariance of the linearly 
divergent term has been checked explicitly in one-loop perturbation theory. The 
same argument can be applied to the evaluation of the subtraction constants 
which appear in the definition of a finite kinetic energy operator, to be discussed 
in the following section. 


It may appear more natural to define A using eqs.(|^) and but with v deter¬ 
mined at a small value of t, {t* say, with \/t* ^ Aqcd)- In particular it may seem 
that the value of A obtained from a measurement at “short distances” can be used 
more reliably to determine some standard short distance mass in QCD say, 

with /r ^ Aqcd) using perturbation theory. This is not the case however, since in ad¬ 
dition to the non-perturbative contribution of 0{Aqcd) to A, there is a perturbative 
contribution which is proportional to 1 /t*, 


i^pert (t ) — 


1 ( |>S'^ert(^.r -Fa;0,0)I 

a 


|5'pert(^. t*; 0 ,o)| 


^ DsCf 
a dvr 
OsCf 7b 
dvr t* 


7b In + ^) + 0{a1) 

+ 0{al) (as a —> 0 ) 


( 2 d) 

(25) 


where 5'pert is the heavy quark propagator in perturbation theory, which for t > 0 takes 
the form 


‘S'pirt(T t; 0) = (^) <5*^ 


1 - (7b ln(t/a) -b Ct) 


(26) 


and the anomalous dimension of the heavy quark field ( 7 .^) and ct are constants (in the 
Landau gauge 7 ^ = — 6 ). Cf is the eigenvalue of the quadratic Casimir operator in the 
fundamental representation (Cir=d/3). In order to determine A, from the propagator 
computed at a finite value of t*, the perturbative contribution must be subtracted. The 
evaluation of the term proportional to 1/t* in perturbation theory (which in practice 
can only be performed up to some low order), becomes less accurate as t* is decreased. 
The reason is that, although the calculation of the coefficient of the term proportional 
to 1/t* becomes more accurate as t* decreases, the presence of the factor 1/t* implies 
that the subtraction becomes larger numerically and that the error due to (unknown) 
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higher order pertubative corrections also increases, reducing the accuracy of the result 


for A. For this reason we propose to define 6m from eq.(21), i.e. from measurements 
of the propagator at large values of t. In some simulations it may not be possible to 
compute the propagators at sufficiently large values of t for a plateau to be reached 
(i.e. for the ratio of the propagators on the right hand side of eq.(|^) to be independent 
of t). In those cases it may be necessary to determine mq from measurements taken 
at intermediate values of t, and to perform the subtraction of the terms proportional 
to 1/t, either by using perturbation theory or by fitting — i/ to a function of t and 
extracting the asymptotic value (i.e. the value as t ^ oo). In refs.|jl5|, |l^ it has been 
shown that the latter method can be used to give a precise determination of A, and 
hence of mq. 


3 The Kinetic Energy - Ai 


In this section we present our proposal for the elimination of power divergences and 
renormalons from the matrix elements of higher dimensional operators. We illustrate 
our method by considering explicitly the kinetic energy operator JiD^h. The matrix 
elements of this operator contain power divergences because it can mix with the lower 
dimensional operators hD^h and hh^. A subtracted kinetic energy operator, one which 
is free of power divergences, is of the form 

hD'gh = hD^h — ci(hDih + 5mhh) — C 2 hh (27) 

where the constants ci and C 2 are fixed by imposing appropriate renormalisation con¬ 
ditions. We propose to define ci and C 2 by imposing that the matrix element of JiD'gh 
between quark states with k = 0 (where k is the momentum of the quark), and in the 
Landau gauge, vanishes 

{h{k = 0)\hDlh\h{k = 0)) = 0. (28) 


Although not unique, this is perhaps the most intuitive definition of the kinetic energy 
of the heavy quark in a hadron. Specifically we determine the constants ci and C 2 by 
using [| 


Cl + C2tx — 


J2x,yJ2ty=0,t^ ( 0 I Hx, tx) h{y, ty)Dlh{y, ty) /i(0,0) | 0) 

4 )^( 0 , 0 ) | 0 ) 


(29) 


® The relation between power divergences in matrix elements of the particular operator hD^h and renor- 
malon ambiguities in dimensional regularisation is subtle and not fully understood. A detailed discussion of 
this subject can be found in ref. Q. 

^From now on we will work in lattice units, setting a = 1. 
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In deriving eq.([^) it is implied that we are using the action C given in eq.(19) which 
contains the residual mass counterterm. However, from the discussion in the previous 


section we can readily see that this is equivalent to using the action C given in eq.(12) 
which has no residual mass term. The only difference in using the actions C and C 
is a factor exp(z/t 3 ;) in both the numerator and denominator of the right hand side of 

eq.®. 

For some important applications it is only the constant C 2 which is required. This 
is because 

-{JiD^h + 6mhh)S'f^{x; y) = 6‘^{x — y ), (30) 


1 


1 + 6m 

so that if the operators in a correlation function are separated (i.e. up to contact 
terms), then the term proportional to ci vanishes. C 2 can also be determined directly 
by eliminating the sum over ty in eq.( 


C 2 = 


J2x,y {0 I Hx, tx) h{y, ty)Dlh{y, ty) /i(0,0) | 0 


y)^y' 


(31) 


J2x{^\hix, tx)H0,0) | 0 ) 

for ty ^ O^tx- 

Having defined the subtracted operator Ai can be determined from a com¬ 

putation of two- and three-point correlation functions in the standard way. Consider 
the meson three-point correlation function (the extension of this discussion to baryons 
is entirely straightforward) 


Cg2^{tx,ty) = ^(0|Jr(f,tx) h{y,ty)Dlh{y,ty)jl{0,0)\0) 

X 

For sufficiently large values of ty and tx — ty 

Cg2^{tx,ty) ^ -■iy)tx) 

where Ai is defined by eq.(§), using the subtracted kinetic energy operator, 

{H\hDlh\H) 


Ai = 


2mH 


(32) 


(33) 


(34) 


and H is the lightest meson state which can be created by the operator Jjt. A convenient 
way to extract Ai is to consider the ratio 


Cg2{tx,ty) 

C{tx) 


Ai 


(35) 


As usual Ai must be evaluated in an interval in which R(tx,ty) is independent of the 
times ty and tx, so that the contribution from the excited states can be neglected. By 
the same argument as was given after eq.(29), the ratio in eq.(35) can be evaluated 


11 










using the action C of eq.(12) with no residual mass term. In the present case the 
difference between the matrix elements of the subtracted and unsubtracted operators 
is given by 

{H\hDlh\H) _ {H\h&h\H) 


Ai = 


C 2 • 


(36) 


2mH 2mH 

We conclude this section by presenting the definition of the subtracted quark mass 
up to, and including the terms of O^l/mq), 

Ai 


Mb = rriQ + S — v — 


2 m| 


(37) 


where Mb is the spin-averaged mass, Mb = \{Mb + ‘?>Mb*)i which has no contribution 
from the chromomagnetic operator. Eq.(|37|) must be modified to include the effects of 
perturbative corrections. We denote the renormalised kinetic energy operator in some 
continuum renormalisation scheme by In one loop perturbation theory we 

have 


^-^cont ^ 


a 


1 + 


from which we derive 


Mb = rug ^ E — V — 


1 + 

47r 


2m^ 


(38) 


(39) 


The term proportional to in eq.(39) is absent in continuum formulations of the 
HQET, and is a manifestation of the lack of reparametrisation invariance in the lattice 
version. It has been calculated in ref. ||l^ . 

Although we have restricted our explicit discussion to the kinetic energy operator, 
clearly the same techniques can be applied to a wide class of operators. This includes, 
for example, the operators whose matrix elements determine the l/mq corrections to 
exclusive leptonic and semileptonic decays of S-mesons ||^. In each case one can con¬ 
struct linear combinations of higher and lower dimensional operators which are free 
of power divergences and renormalon ambiguities, by imposing appropriate normali¬ 
sation conditions for matrix elements between quark states in a fixed gauge, and at 
given momenta. This approach is a particular application of the general method for 
the non-perturbative normalisation of lattice operators proposed in ref.||lf 


4 Matching 

In the proceeding sections we have proposed a method for defining higher-dimensional 
operators whose matrix elements are free of (ultra-violet) renormalon ambigui¬ 

ties and power divergences, and which can be computed in lattice simulations. In order 
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to derive physical predictions from these matrix elements it is necessary to calculate 
the corresponding coefficient functions, i.e. the Cn,a of eq.(j^. The subtraction of the 
ultra-violet renormalons from the matrix elements, implies the elimination of the cor¬ 
responding infra-red renormalons from the coefficient functions, which can therefore be 


computed in perturbation theory [17, M, 21|. To illustrate the “matching” procedure 


consider a simple situation for which we can write eq.(|^) as: 


0«™(M) = C’i(^, + -hc2(^. 4)02''®'"”(A) + 0(l/m|) (40) 


A 


M ’ M 




M ’ M' 


i.e. where there is a single operator in each of the first two terms of the heavy quark 
expansion (the discussion can readily be extended to the general case). In eq.(p0[) we 
have exhibited explicitly the dependence on the scale M used to normalise the QCD 
operators (i.e. the operators in the full theory). A is the (hard) ultra-violet 

cut-off in the effective theory (for example in the lattice formulation of the HQET 
A = a“^), and the operators on the right-hand side of eq.(^) are bare operators, uiq 
can be considered as the subtracted pole mass mg defined in section 

Note that although we have chosen to write eq.(^0[) in terms of the bare operators 
in the effective theory, this equation is equivalent to eq. i), which is written using 
renormalised operators. From the bare operators (computed in lattice simulations for 
example), one can determine the corresponding ones renormalised in a given (contin¬ 
uum) scheme by using perturbation theory, or by some non-perturbative method (such 
as that given in ref.[l^). The determination of the matrix elements of renormalised 
operators of the effective theory may be a convenient intermediate step, but it is not 
necessary. We choose instead to compute the matrix elements of the QCD operators 
qQCD from the bare operators of the effective theory. 

With a hard cut-off the Borel transforms of the coefficient function Ci and the 
matrix elements of do not have renormalon singularities at u = 1/2. For 

example, in the large Nf limit, the Borel transform of Ci in the vicinity of u = 1/2 
will have the structure 


- mg A 1 



A /A\-2“ 
mg \m) 


(41) 


where the tilde denotes the Borel transform. Although the residue of the pole vanishes 
at u = 1/2, both Ci and the matrix elements of contain terms which, in 

perturbation theory, diverge linearly with the ultra-violet cut-off A. The use of a 
hard factorisation scale to organise the heavy quark expansion has been suggested by 
Bigi et al. § (see also ref. j^ for a very recent study of the definition of the higher 
twist operators relevant for studies of deep inelastic structure functions). Our proposal 
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is to subtract these power divergences non-perturbatively, and to use the subtracted 
operators, which are free of power divergences and renormalon ambiguities, as the basis 
for the expansion in eq.(|I|). 

Following the discussion in the preceding sections, we define a subtracted operator 
which does not mix with 

(42) 

WT' / 

where the dependence of the matrix elements of 02-5 (^) ^ most logarithmic. 

c(A) is computed non-perturbatively, whereas the functions Ci and C 2 are calculated 
in perturbation theory. Eq.(40) can be rewritten as: 

A 






^ ,mQ A c{K) niQ 

M" mg M ’ M‘ 


O 


HQET 
1 


(A) 


1 

mg 


+ :f7Cj(^,A)0»«'=’-(A)+0(l/m|) (43) 


By using a combination of perturbative and non-perturbative techniques, we have en¬ 
sured that the coefficient of in eq.(43) does not contain terms which diverge 

linearly with A, nor ambiguities due to renormalon singularities at u = 1/2. The 
remaining ambiguities are of 0{AQ^j^/mQ) or less, and are associated with the renor- 
malons at n = 1, | • • •, or the corresponding power divergences in matrix elements of 
operators of higher dimension than These can be eliminated by generalising 

our procedure to higher orders of the heavy quark expansion, by defining subtracted 


operators 


O 


HQET 

4;S 


which cannot mix with lower dimensional ones, and 


using these subtracted operators as the basis of the expansion in eq.(||). Indeed by 
using such a basis one eliminates the ambiguities from all coefficient functions, up 
to the order in Ijmg for which the subtraction coefficients have been computed. Of 
course in general the coefficient function of depends logarithmically on A, the 

dependence being given by the anomalous dimension of 

Throughout this paper we have been stressing the necessity of performing the sub¬ 
tractions of power divergences non-perturbatively. We now present a specific example 
demonstrating this explicitly. Consider the evaluation of the subtraction constant c(A) 
in the case for which = Jih and = hv ■ D h as happens in the expansion 

of the QCD operator QQ, where Q is the field of the heavy quark in the full theory, and 


in the discussion of A in section 2.2. The Borel transform of the forward matrix element 


of between heavy quark states of momentum k {T = {h{k) \ hv ■ Dh \h{k))), i 

the vicinity of tt = 1/2, has the structure 


m 


T(n) (X 


1 


(1 - 2u) 


V ■ k 


V ■ k 

~A' 


—2u 


- A 


(44) 
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The residue of the pole vanishes at u = 1/2, in a similar way to that in the coefficient 
function Ci in eq.(41). The two terms in the square brackets in eq.(44), contribute 
equal and opposite ambiguous terms of 0{Aqcd) to the matrix element of at 

V ■ k = 0. In perturbation theory however, the first term does not appear at n • A; = 0, 
and hence the Borel transform of c(A) contains a renormalon at n = 1/2. It is for 
this reason that we insist on the non-perturbative determination of the subtraction 
coefficients. 

One can also imagine performing the matching completely non-perturbatively, by 
simulating both the HQET and QCD (with a bottom quark) on the lattice. However 
the latter requires a very small lattice spacing, a~^ S> mh, which will not be possible 
for some time to come. Moreover, once one is able to simulate the 6-quark directly 
(and reliably) on the lattice, the necessity of using the HQET is removed. It may still 
however, be a useful guide to scaling properties and symmetry relations. 


5 Conclusions 

In this paper we have presented a method for defining higher dimensional operators 
of the HQET, in such a way that their matrix elements are free of ambiguities due to 
(ultra-violet) renormalon singularities, and of power divergences. We have illustrated 
our approach by proposing a physical definition of A (= mn — w-g) and of the matrix 
elements of the subtracted kinetic energy operator {{H\Dg\H)). The dehnition of the 
higher dimensional operators involves the snbtraction of lower dimensional ones with 
the same quantum numbers. The subtraction coefficients are determined by imposing 
normalisation conditions on Green functions between quark and gluon states (or in the 
case of A, on the heavy quark propagator). Lattice simulations of the HQET allow for a 
numerical evaluation of the snbtraction coefficients, as well as of the matrix elements of 
the subtracted operators. The renormalisation procedure proposed above can, however, 
be applied with any other non-perturbative method for computing matrix elements in 
effective theories. Our approach can also be extended to other cases, for example to 
the higher twist contributions to the structure functions of deep inelastic scattering. 

Having defined operators 0^2^"^ whose matrix elements are free of renormalon 
ambiguities, we still have to match the HQET operators onto those of QCD, i.e. to 
determine the coefficient functions Cn,Q of eq.(|^. This is done nsing perturbation 
theory, by calculating the matrix elements of and between suitable 

external states, and combining the results with the subtraction coefficients which have 
been computed non-perturbatively (see section above). 
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The procedures defined in this paper allow one to study quantitatively many im¬ 
portant physical processes and quantities in heavy quark physics using a systematic 
expansion in the mass of the heavy quark. These include the leptonic and semi-leptonic 
decays of heavy mesons and baryons, as well as relations between the masses and life¬ 
times of heavy hadrons. The “subtracted pole mass” is a suitable parameter for the 
expansion. 
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